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An  asymptotic  theory  is  presented  to  analyze  perturba*^ 
tions  of  bifurcations  of  the  solutions  of  nonlinear  problems. 
The  perturbations  may  result  from  imperfections,  impurities, 
or  other  inhomogeneities  in  the  corresponding  physical 
problem.  It  is  shown  that  for  a wide  class  of  problems  the 
perturbations  are  singular.  The  method  of  matched  asymptotic 
expansions  is  used  to  obtain  asymptotic  expansions  of  the 
solutions.  Global  representations  of  the  solutions  of  the 
perturbed  problem  are  obtained  when  the  bifurcation  solutions 
are  known  globally.  This  procedure  also  gives  a quantitative 
method  for  analyzing  singularities  of  nonlinear  mappings  and 
their  unfoldings.  Applications  are  given  to  a simple 
elasticity  problem,  and  to  nonlinear  boundary  value  problems. 
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1.  Introduction. 

Bifurcation  theory  is  a study  of  the  branching  of  solu- 
tions of  nonlinear  equations, 

(1.1)  fly; A)  “ 0 , 

where  f is  a nonlinear  operator,  y is  the  solution  vector, 
and  A is  a parameter.  It  is  of  particular  interest  in  bifur- 
cation theory  to  study  how  the  solutions  y(A)  of  (1.1)  and 
their  multiplicity  change  as  A varies.  Thus  we  refer  to  A as 
the  bifurcation  parameter.  A bifurcation  point 

of  a solution  branch  y(A)  of  (1.1)  is  a point 
l XQ # y < Aq) 3 from  which  another  solution  y^(A)  branches.  That 
is,  y(A©)  “ y^A©) , and  y(A)  ? y^(\)  for  all  A in  an  interval 
about  A©.  Thus  at  a bifurcation  point  there  is  a tran- 

sition in  solution  multiplicity. 

Typical  applications  of  bifurcation  theory  are  in  elastic 
and  hydrodynamic  stability.  In  elastic  stability  theory,  bi- 
furcation is  called  buckling,  A represents  an  applied  load, 
and  the  bifurcation  points  are  called  buckling  loads.  In  hy- 
drodynamic stability,  bifurcations  are  called  transitions,  A 
is  a flow  parameter  such  as  the  Reynolds  number,  and  bifurca- 
tion points  are  called  transition  or  critical  points. 

In  experiments  and  in  real  applications,  the  sharp  tran- 
sitions of  bifurcation  rarely  occur.  Small  imperfections,  im- 
purities or  other  inhomogeneities  tend  to  distort  these  tran- 
sitions. For  example,  in  the  buckling  of  thin  elastic  rods, 


3 


plates,  and  shells,  the  classical  mathematical  theory  usually 
results  in  a bifurcation  problem.  However,  experiments  show  that 
the  transitions  are  smooth  rather  than  sharp.  The  discrepancy  is 
attributed  to  small  deviations  in  the  shapes  of  the  experimental 
specimens  from  the  perfect  shapes  assumed  in  the  theory,  and 
to  other  physical  and  experimental  imperfections.  Indeed,  it 
has  been  demonstrated  experimentally  that  when  the  magnitudes 
of  the  imperfections  are  decreased,  the  experimental  results 
approach  the  predictions  of  the  bifurcation  problem.  Moreover, 
it  has  been  demonstrated  by  experiments  and  by  approximate 
calculations  on  "imperfection-sensitive"  elastic  structures 
that  small  imperfections  may  have  a large  effect. 

Small  experimental  imperfections  and  inhomogeneities  are 
of  significance  in  other  nonlinear  problems,  such  as  hydro- 
dynamic  stability,  as  we  shall  show  in  a future  publication. 

To  analyze  the  perturbation  of  bifurcations  produced  by 
small  inhomogeneities,  we  shall  consider  the  nonlinear  problem 

(1.2)  F[y;X,6]  - 0 , 

where  the  additional  parameter 

< characterizes  the  magnitude  of  the  inhomogeneities.  We 
assume  that  as  £ ♦ 0,  (1.2)  and  its  solutions  y(X,6)  re- 
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y(X,0)  - y (X)  , 

(1.3) 

F[y(X,0)jX,0]  •=  f ly  (X)  j X]  . 

We  refer  to  (1.2)  as  the  perturbed  problem,  and  (1.1)  as  the 
bifurcation  problem.  Furthermore,  we  refer  to  the  solution 
branches  of  (1.2)  and  (1.1)  as  the  perturbed  and  bifurcation 
branches , respectively. 

In  this  paper,  the  method  of  matched  asymp- 

totic expansions,  see  e.g.  [1,2,3],  is  used  to  obtain  formal 
asymptotic  expansions  of  the  solutions  of  the  perturbed  prob- 
lem (1.2)  as  6 -*>  0.  To  motivate  our  analysis,  we  consider  the 
simple  algebraic  equation, 

(1.4)  F [y ; X , 6 ] = y2  + (l-X)y  - X6  « 0 . 

The  solutions  of  the  corresponding  bifurcation  problem 
(6  - 0)  are, 

(1.5)  y - Yq(X)  = 0 , y - Y^X)  = X - 1 . 

They  are  sketched  in  Figure  lb.  Thus,  X ■ 1,  y - 0 is  the 
only  bifurcation  point  of  the  solutions  of  the  bifurcation 
problem.  The  solution  branches  of  the  perturbed  problem  (1.4) 

are 

y - y±(x,«)  5 ^ ± J[(X-d2+4X«i1/2  . 


(1.6) 
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They  are  sketched  in  Figures  la  and  lc  for  two  representative 
small  positive  and  small  negative  values  of  6,  respectively. 

We  observe  that  as  j + 0,  the  perturbed  branches  (1.6)  approach 
the  bifurcation  branches  (1.5).  Since  the  perturbed  branches 
are  smooth  and  the  bifurcation  branches  have  a "corner"  at 
the  bifurcation  point,  the  approach  as  6 0 is  non-uniform 

in  any  X interval  containing  the  bifurcation  point.  Thus  (1.4) 
is  a singular  perturbation  of  the  bifurcation  problem.  This 
suggests  the  use  of  the  method  of  matched  asymptotic  expansions. 

The  outer  expansions  of  the  perturbed  branches  are  regu- 
lar perturbations  in  6 of  the  bifurcation  branches.  As  we 
shall  see,  these  expansions  are  not  valid  in  neighborhoods  of 
bifurcation  points  or  other  singular  points  of  the  bifurca- 
tion branches.  (Singular  points  are  defined  in  Section  2.) 

In  these  neighborhoods,  inner,  or  boundary  layer,  expansions 
are  used  to  represent  the  perturbed  branches.  A matching 
procedure  is  then  used  to  "connect"  the  inner  and  outer  expan- 
sions. Uniform  asymptotic  expansions  of  the  perturbed 
branches  can  then  be  obtained  by  forming  a composite  expansion; 
see,  e.g.  [1,2,3], 
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The  asymptotic  expansions  that  we  obtain  are  valid  for 
any  ranges  of  X for  which  the  bifurcation  branches  are  known, 
either  exactly,  or  approximately,  for  example,  by  other  per- 
turbation expansions  or  by  numerical  computations.  Thus,  when 
bifurcation  branches  are  known  globally,  our  method  gives  global 
repres jntations  of  the  solutions  of  the  perturbed  problem.  Pre- 
vious studies  [4-8]  of  such  perturbations  of  bifurcations  are 
local  analyses  valid  only  near  bifurcation  points. 

In  Section  2,  we  apply  the  method  to  a general  nonlinear 
equation  (1.2).  Relatively  mild  restrictions  are  imposed  on 
the  operator  F.  We  show  how  the  determination  of  the  inner 
expansions  is  reduced  to  the  analysis  of  algebraic  eqv  ons. 
This  is  not  surprising  since  the  Lyapunov-Schmidt  [9] , and 
other  theories,  effect  a similar  reduction  in  classical  bi- 
furcation theory.  We  emphasize  that  the  analysis  of  these 
algebraic  equations  is  generic  to  a general  class  of  nonlinear 
equations.  This  is  illustrated  in  Section  4,  where  we  apply 
the  method  to  nonlinear  boundary  value  problems. 
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2.  A General  Theory. 

We  assume  that,  the  solutions  y(  X , 6 ) of  the  perturbed 
problem  (1.2)  are  elements  of  a real  Hilbert  space,  F is  a 
nonlinear  operator  on  that  space,  and  X and  6 are  real  valued 
parameters.  For  simplicity  of  notation,  the  possible  depen- 
dence of  F and  y on  other  variables  is  suppressed.  We  shall  use 
the  bracket  notation  <*,*>  for  the  inner  product.  Sufficient 
conditions  on  F for  the  application  of  our  method,  other  than 
obvious  differentiability  requirements,  will  be  stated  as 
needed. 

The  bifurcation  problem  corresponding  to  (1.2)  is 

(2.1)  F[y (X,0) ?X,0]  = f [y (A) , A]  » 0 . 

We  assume  that  the  bifurcation  points  of  a solution  y(A)  of 

(2.1)  are  determined  from  the  linear  problem, 

(2.2)  fyly(A)  ;A]4>  - 0 

that  is  obtained  by  linearizing  (2.1)  about  y ( A) . The  linear 
operator  fy(y(A);A]  is  the  Frechet  derivative  of  f evaluated 
at  y(A)  and  X.  The  values  A ■ Ap  for  which  (2.2)  has  solu- 
tions ♦ jf  0,  and  the  corresponding  values  y(XQ),  define  the 
singular  points  of  the  solution  branch  y(A).  In  particular, 

bifurcation  points  are  singular  points.  For  many  problems 
all  the  singular  points  are  bifurcation  points. 
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Por  the  analysis  that  follows,  we  define  the  sub-branches 
of  a bifurcation  branch  y(l)  as  the  segments  of  y(l)  connecting 
two  successive  singular  points  of  the  branch.  In  addition,  we 
shall  consider  the  segments  of  a bifurcation  branch  that  connect 
the  "last"  singular  points  with  the  end  points  of  the  branch, 
which  may  be  at  *»,  as  sub-branches . 


The  Outer  Expansions. 


Let  y = y^ (X)  denote  a sub-branch.  We 
seek  asymptotic  expansions  as  6 0 of  the  solutions  y(X,6) 

of  the  perturbed  problem  (1.2)  near  y^fX)  in  the  form, 


(2.3) 


y(x,«) 


yjU)6 


j 


The  coefficients  y^,y2*...#  are  determined  by  inserting  (2.3) 
into  (1.2)  and  equating  to  zero  the  resulting  coefficient  of 
each  power  of  6.  This  leads  to  the  following  recursive  sys- 
tem of  equations 


(2.4)  F®yx (X)  - -F®  £ Rjd) 

(2.5)  pJy2(X)  - - £ [(P^yj.)?!  + ZFyjY!  + p|51  = Rj (» 


and  in  general 

Fj y^X)  - Rj  (X) 


(2.6) 
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The  derivatives  of  F with  respect  to  y are  functional,  or 

Frechet,  derivatives.  The  superscript  zero  denotes  the  corresponding 

operator  evaluated  at  6 * 0,  e.g.  F°  = F°(X)  = py  Iy0 (^) ; • 

The  inhomogeneous  terms  R^  depend  on  y^y^,  • • • 'yj-i* 

The  solutions  of  (2.4)-(2.6)  are  written  as 

(2.7)  Yj  = (Fy)’1^  * j - 1.2 

where  (F°)_1  is  the  inverse  of  F°.  The  solutions  (2.7)  exist 
for  all  X that  are  not  singular  points  of  Yq(X).  At  the  sin- 
gular points,  F°(X)  = fy[yQ(X);X]  is  not  invertible.  The  co- 
efficients Yj  exist  at  a singular  point  XQ  of  yQ(X),  only  if 
the  R..  (XQ)  are  orthogonal  to  the  null  space  of  the  adjoint 
operator  (F°(XC))*.  We  assume  that  the  operator  F is  such 
that  this  condition  is  violated  for  some  j.  In  particular  we 
take  j = 1.  That  is,  we  assume  that  there  is  an  element  4>* 
of  the  null  space  of  F°*  such  that 

(2.8)  I = <R1(X0),<^*>  = -<Fg(X0),$*>  ? 0 . 

Then,  as  we  shall  show,  c.f.  Figure  1,  (2.8)  implies  that  the  in- 
homogeneities "destroy"  the  bifurcations.  It  is  this  type  of 
perturbation  of  bifurcation  that  we  shall  study  in  this  paper. 

The  condition  (2.8)  is  satisfied  for  the  problems  analyzed  in 
Sections  3 and  4,  and  for  many  other  problems. 

If  I « 0,  a similar  analysis,  which  we  do  not  present  in 
this  paper,  is  required.  Then  the  inhomogeneities  need  not 
destroy  bifurcation;  the  perturbed  problem  may  be  a bifurcation 
problem.  An  example  of  this  is  discussed  in  [10).  See  [11,12] 
for  other  perturbations  of  bifurcations  that  result  in  secondary 
bifurcation. 
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The  expansion  (2.3),  which  is  not  valid  at  the  singular 
points  of  Yq(M»  is  called  the  outer  expansion  corresponding 
to  Yq(M.  The  Yj(X),  j « 1,2,...,  are  called  the  outer 
coefficients. 


The  Inner  Expansions. 

We  obtain  expansions  that  are  valid  near  \Q  as  6 -*■  0, 
by  "stretching"  the  neighborhoods  of  XQ  through  the  transfor- 
mation 

m 

(2.9)  X « XQ  + £ya  + ^ ^(V*)1  , 

where  a > 0 is  a constant,  and  the  small  parameter 

)j  is  defined  by 

(2.10)  fi(y)  - (sgn6) yb  . 

The  constant  b in  (2.10)  must  be  positive  because 

we  require  6 (0)  « 0.  The  values  of  a and  b are  de- 

termined by  the  nonlinearity  of  F near  XQ,  as  we  shall  demon- 
strate. 

Denoting  the  dependent  variable  for  the  analysis  of  the 

inner  expansion  by  z(y)  «*  y(X(y)  ,6  (y ) ) , we  seek  asymptotic 

expansions  of  the  solutions  of  the  perturbed  problem  (1.2) 

near  XA  in  the  form 
o 
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(2.11) 


ZjM^  f 


where  X = XQ,  zq  “ Yq (*q)  *s  ^e  ®in9ular  point.  The  expan- 
sions (2.11)  are  called  the  inner  expansions, 

and  the  z ..  are  called  the  inner  coefficients.  We  insert 
(2.9)- (2.11)  into  (1.2)  and  equate  to  zero  the  coefficient  of 
each  power  of  y.  This  leads  to  the  following  system  of  equa- 
tions to  successively  determine  the  zy. 

-lO_  _ —O.O  nO/1  » iO  — 


(2.12)  F^  - - Fj (*„>«“  = Pj.  . 

(2.13)  F°z2  - - |<(F”2z1)z1  + + (X°)2F°X 


+ 02FzS*l  + 22uFX«  + 2 p2  - 


and  in  general, 

(2.14)  F°Zj  - Pj  . 

Subscripts  on  X and  6 denote  derivatives  with  respect  to  y. 
The  superscript  zero  means  that  the  quantities  are  evaluated 
at  y ■ 0.  Since  X(0)  - XQ  and  6(0)  ■»  0,  the  operators  in 
(2.12)- (2.14)  are  evaluated  at  the  singular  point.  Thus,  for 
example, 

- F2tz(0);Xo,0)  - F°(Xo)  , 

(2.15) 

Fj(X0)  - F6(*(0),Xo,0J  . 


I 
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We  require  that  the  derivatives  of  X and  6,  evaluated  at 
y » 0,  are  bounded  to  insure  that  the  right  sides  of  (2.12)- 
(2.14)  are  bounded.  Then  (2.9)  and  (2.10)  imply  that  a and  b 
are  positive  integers.  Furthermore,  we  observe  that  if  b *=  1, 

then  the  inner  expansion  (2.11)  is  a power  series  in  6. 

* 

We  shall  not  consider  such  expansions  in  this  paper.  Thus 
b 5^  2,  and  hence, 

(2.16)  - 0 . 

Since  F®  is  the  linearized  operator  at  the  singular 
point,  the  inhomogeneous  equations  (2.12)- (2.14)  have  solu- 
tions if  and  only  if  each  is  orthogonal  to  the  null  space 

of  F°*.  We  assume  that  the  null  space  is  one-dimensional,  so 
z 

that 

(2.17)  "0»  j ■ 1,2,...  . 

We  shall  consider  multi-dimensional  null  spaces  in  a subse- 
quent publication.  The  solvability  condition  (2.17)  with 
j ■ 1,  and  equation  (2.16)  imply  that 

(2.18)  *l<*l'**>  - 0 • 

Thus  if  the  projection 

(2.19)  A £ <*£,♦*>  * 0 , 

then  (2.18)  implies  that  X®  ■ 0,  and  (2.12)  implies  that 

"■quation  (2.17)  with  j - 1,  and  assumption  (2.8)  imply 
that  the  expansion  (2.9)  for  X is  a power  series  in  6. 
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(2.20)  z1  « A*  , 

where  the  eigenfunction  $ is  normalized  by  “ 1*  The 

amplitude  A is  to  be  determined.  If  A * 0,  then 

(2.21)  z1  - A*  + X°ZX  , where  Zx  = -(F°)_1F°  . 
Quadratic  Nonlinearities. 

Condition  (2.17)  with  j - 2,  (2.16),  and  (2. 19) -(2. 21) 
imply  that  A must  satisfy  the  quadratic  equation 

(2.22a)  <(F°Z*)*,**>A2 

+ + * 0 • 

if  A ^ 0,  and 

(2.22b)  <(F°Z*)*,**>  A2  + <U»J1*)*1+(fJ1*1)*+2»Jx»),*»»xJx 

+ < < x°>  2<  t (f^)  V2F*x  Vr!  x]  ,♦•> 

+ * 0 ' 

if  A * 0.  Equations  (2.22)  are  the  analogues  of  the  bifurca- 
tion equation  of  standard  bifurcation  theory.  The  solutions 
of  these  equations  determine  A as  a function  of  X and  of  the 
inhomogeneity,  for  X near  XQ.  For  each  real  root,  we  obtain 
a coefficient  z^t  which  gives  the  leading  term  in  an  inner 
expansion. 
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If 

(2.23) 


Q = <(?„♦)♦»♦*>  + 0 , 


then  in  order  for  the  solutions  A of  (2.22)  to  depend  on  both 
the  bifurcation  and  inhomogeneity  parameters,  we  must  have 

if  A + 0 , 


(2.24)  + 0 , and  \ 


XZ»  " 0 • 


* 0 - 


if  A - 0 . 


Since  a and  b are  positive  integers,  this  implies  that 
6 » (sgn6)y2  (b  » 2),  and  X ■ + £y2  + 0(y3)  (a  - 2)  if 

A / 0,  and  X - XQ  + £y  + 0(y2)  (a  = 1)  if  A « 0.  This  is 
equivalent  to 


(2.25a) 


X 


+ €|6|  + 0 (62)  , for  AM; 
X©  + ?|6|1/2  + 0(6),  for  A « 0 ? 


(2.25b)  s 


*0  + A(C)*|6|1/2  + 0(6)  , for  A ? 0 , 

*0  ♦ (A(C)4+CZ1)|6|1/2  + 0(6)  , for  A - 0 . 


Equations  (2.25),  with  the  amplitude  A determined  by 
(2.22),  give  parametric  descriptions,  with  £ as  the  parameter, 
of  the  solutions  of  the  perturbed  problem  near  the  singular 
point.  Condition  (2.23)  characterizes  the  nonlinearity  near 
XQ  as  quadratic.  The  two  types  of  quadratic  nonlinearity 

are  in  turn  characterised  by  A » 0 or  A j*  0. 


% 


Cl 
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Higher  order  terms  in  the  expansions  can  be  obtained  by 
analyzing  (2.13)  and  (2.14).  The  amplitudes  Aj  of  the  eigen** 
function  $ in  each  Zj  for  j > 1,  are  then  solutions  of  linear 
inhomogeneous  algebraic  equations.  The  solvability  condi- 
tions for  these  equations  determine  the  parameters  for 
j > 2. 


Cubic  Nonlinearities. 

If  Q - 0,  then  Equations  (2.22)  imply  that 

- 0 , if  A^O 

1°  - 0 , if  A - 0 . 

Otherwise,  for  (2.22a),  the  variable  £ ■ 1/2  Ayy  would  be  a 

fixed  number,  and  for  (2.22b)  there  would  be  no  A which  is 
bounded  as  A **■  Aq,  and  simultaneously  depends  on  the  inhomo- 
geneity. Hence  Equations  (2.22)  are  satisfied  identically 
and  (2.13)  gives 

(2.27)  z2  - Aji  + A222  ♦ A°yZ2  , 

where  A2  is  to  be  determined,  and  Z2  and  Z2  are  defined  by 

(2.28)  *2  * " # *2  5 “ I(P^)1fA  * 

The  solvability  condition  (2.17)  with  j ■ 3 implies 

that. 


(2.26)  6°  - 0 , and  ^ 


I 


(2.29a) 


CA3  + (X°  - 0 » 


UVV  X‘ 


if  A f 0 , 


(2.29b)  CA3  + <{2F®zZ2«+4(F°2<|.)Z2+3F°x^)^*>X°i|A 


+ 6°  <F?,4>*>  *=  0 , 

yyy  6 


if  A * 0 , 


where  the  coefficient  C is  defined  by. 


(2.30)  C = <<((F®ZZ4>)*)*  + 4(F°z*)Z2  + 2(FzzZ2)^)'^*>  * 


These  are  cubic  algebraic  equations  for  A if  C f 0.  Notice 
that  the  A2  term  is  missing  in  both  equations.  For  each 
real  root  of  (2.29)  we  obtain  an  inner  expansion.  Thus 
there  is  always  at  least  one  inner  expansion.  The  condition 
C j*  0 characterizes  cubic  nonlinearities. 

If  C f 0,  then,  as  in  the  analysis  of  (2.22),  in  order 
for  A to  depend  on  both  the  bifurcation  and  inhomogeneity 
parameters,  v/e  must  have 


(2.31)  6°w  * 0 , and 


Sum  * 0 ' if  A * 0 ' 


Su  * 0 ' lf  * " 0 • 


Then  { ■ (sgn6)y3,  X - ♦ $y3  + 0(y*)  if  A + 0 and 

X ■ 1 ♦ Ey2  ♦ 0(y3)  if  A - 0,  or  equivalently 
o 
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x - xQ  + *1*1  ♦ °<*4/3> 

(2.32)  l 

z - z0  + A«  | £ | 1/3  + 0(62/3) 


for  A i*  0 , 


X » XQ  + 5|£|2/3  + 0(6) 

(2.33) 

z - zQ  + A+ | 6 | 1/3  + 0(62/3) 


for  A “ 0 • 


If  C = 0« we  can  show  that  the  solvability  condition 
<P3,<|>*>  **  0 is  satisfied  identically.  Then  we  must  consider 
solvability  conditions  (2.17)  with  j > 3 to  evaluate  A.  In 

this  way,  we  can  determine  the  values  of  a and  b,  and  we  can 
formally  classify  the  singularities  of  a general  nonlinear 
equation.  This  is  a by-product  of  our  asymptotic  procedure 
for  the  construction  of  solutions  of  the  perturbed  problem. 
Thom  [8]  has  rigorously  classified  the  singularities  of  a 
special  class  of  equations,  using  different  methods. 


We  observe  from  (2.25a)  and  (2.33)  that  for  A * 0, 

X - Xc  ■ 0 (A1^2)  for  the  quadratic  nonlinearity  and 
X - X ■ 0(62^3)  for  the  cubic  nonlinearity.  Thus  the  qua- 
dratic nonlinearity  is  a more  "severe”  perturbation. 


Matching. 

We  have  obtained  outer  expansions  (2.3)  of  the  solutions 
of  the  perturbed  problem  (1.2)  corresponding  to  each  solution 
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problem 

sub-branch  Yq(X)  of  the  bifurcation/ (2. 2) . These  expansions 
are  valid  away  from  the  singular  points  of  y^.  In  addition, 
one  or  more  inner  expansions  (2.11)  have  been  obtained  near 
each  singular  point.  It  is  necessary  to  "connect"  each  of 
the  outer  expansions  to  the  appropriate  inner  expansions,  to 
obtain  uniform  representations  of  the  solutions  of  the  per- 
turbed problem.  For  many  problems  these  connections  are  ob- 
vious, once  the  inner  and  outer  expansions  have  been  obtained 
explicitly.  These  connections  can  be  obtained  systematically 
by  appealing  to  the  matching  procedure  of  the  method  of 
matched  asymptotic  expansions. 


We  derive  the  matching  conditions  by  first  considering 
the  outer  coefficients  near  XQ,  y^(X)  * y^  (X  (y) ),  where 

X(v)  is  given  by  the  stretching  (2.9).  Then  the  inner  expan- 
sion (2.3)  becomes 

m 

(2.34)  y - yj(X(y))  (6gn«)Vb  , 


where  we  have  used  4 « (sgn6)yb.  We  assume  that  for  X near 
XQ,  the  sum  (2.34)  can  be  expanded  as 
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Thus  (2.35)  is  the  outer  expansion  near  XQ.  It  is  expressed 
in  terms  of  the  "inner  variables".  The  outer  coefficients 
depend  on  g,^****  • 

The  outer  expansions  (2.3),  or  equivalently  (2.35),  are 
valid  for  X bounded  away  from  XQ.  The  inner  expansions 
(2.11)  are  valid  for  X near  XQ.  We  assume  that  there  is  a 
common  interval  in  which  both  the  inner  and  outer  expansions 
are  valid;  and  furthermore,  this  interval  shrinks  in  width 
and  approaches  XQ  as  6 0.  Since  the  inner  and  outer  expan- 

sions are  asymptotic  expansions  of  the  same  function,  their 
difference  is  asymptotic  to  zero  in  the  common  interval . By 
using  (2.11)  and  (2.35)  this  gives 


(2.36) 


(nlc-2]c)y  - o 


in  the  common  interval.  Then  it  follows  from,  (2.36),  the 
definition  of  an  asymptotic  expansion,  and  from  (2.9),  that 


(2.37)  lim  tnv-*v]  - 0 , 

K I4" 


for  k « 1,2, 


We  refer  to  (2.37)  as  the  matching  conditions. 

Composite  expansions  can  be  formed  from  the  inner  and 
outer  expansions,  as  is  customary  in  the  method  of  matched 
asymptotic  expansions,  to  give  asymptotic  expansions  in  6 
that  are  uniform  in  X.  We  shall  not  present  these  composite 
expansions. 


J 
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3.  An  Algebraic  Problem. 

As  a simple  illustration  of  our  method,  we  consider  the 
algebraic  equation 

(3.1)  F[y;X,6]  = y3  + 20y2  + (l-X)y  - X6  = 0 , 

where  3 > 0 is  a given  constant.  The  solutions  of  this  equa- 
tion give  the  equilibrium  positions  of  a one-degree-of- 
freedom  model  for  the  buckling  of  an  elastic  arch,  as  we  show 
in  Appendix  I.  The  bifurcation  parameter  X,  which  is  propor- 
tional to  the  load  on  the  arch,  the  inhomogeneity  parameter 
6,  which  is  proportional  to  the  initial  imperfection  in  the 
arch,  and  the  constant  6 are  defined  in  Equation  (A. 2)  of  the 
Appendix.  The  solutions  of  (3.1)  are  easily  obtained,  thus 
providing  a check  of  our  method. 

The  bifurcation  problem  (6  *=  0)  corresponding  to  (2.1) 
is 

(3.2)  P[y;X,0)  - f ty,X]  - y3  + 20y2  + (l-X)y  = 0 . 

The  solutions  of  (3.2)  are 

(3.3)  y - Yq(X)  = 0 , y - Yj(X)  - -P  ± (X-XJl)1/2  , 

They  are  sketched  in  Figure  2b.  The  singular  points  of  these 
bifurcation  branches  are  determined  from  the  linearized  prob- 
lem (2.2). 

Thus  for  y - YQ(X)  = 0,  the  only  singular  point  is  at  X ■ 1, 
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y ■ 0.  It  is  a bifurcation  point,  as  we  see  from  Figure  2b. 
From  (2.2),  (3.2)  and  (3.3),  we  deduce  that  the  singular 
points  of  Y^  are,  the  bifurcation  point  at  X = 1,  and  the 
■limit"  point  at  X *=  X£,  y “ -0.  Furthermore,  this  limit 
point  is  the  only  singular  point  of  Yj, . The  sub-branches 
of  Yc  are  denoted  by  Y*,  and  the  sub-branches  of  Y£  by  Y*+ 
and  Y^~,  as  shown  in  Figure  2b. 

Corresponding  to  each  of  these  five  sub-branches  there 
are  outer  expansions  (2.3)  of  the  solutions  of  the  perturbed 
problem  (3.1).  In  an  obvious  notation,  and  for  the  appro- 
priate range  of  X values,  they  are  given  by 


(3.4) 


y*0)  " 6 + 0(62)  ' 

yU)  “ y(i)  “ Y\l)  + (V2){X-X£-0(X-X£)1/2}"16 

+ 0 (62)  , 

yU)  “ Yi  + (V2)(X-X£+B(X-X£)1/2}"16  + 0(62)  . 


X XX  J,— 

He  observe  that  Y(o)»  y(l)  and  y(l)  are  un^oun<3e^  at  the  bi- 
furcation point,  and  y^j  and  y"^  are  unbounded  at  the  limit 
point.  They  are  shown  by  the  solid  curves  in  Figures  2a  and 
2c.  Since  F^  » -X,  and  the  normalized  eigenfunctions  for 
both  singular  points  are  4 ■ 1,  we  conclude  that  condition 
(2.8)  is  satisfied  at  both  singular  points. 

To  determine  the  nature  of  the  nonlinearity  at  the 
singular  points,  we  evaluate  Q in  (2.23).  Since  F>z  ■ 6z  + 4$, 


t 


we  see  that  Q / 0 at  the  bifurcation  (z  *=  0)  and  limit 
(z  « -0)  points.  Consequently,  the  perturbed  problem  has 
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quadratic  nonlinearities  at  both  singular  points.  Fur- 

thermore, *«  -y,  so  that  A * 0 at  the  bifurcation  point, 
and  A ■ 0 ? 0 at  the  limit  point. 

Thus  at  the  bifurcation  point  XQ  = 1,  y ■ 0,  we  have 
a » 1,  b = 2 and  from  (2.25), 

(3.5a)  A * 1 + £vi  + 0(y2)  « 1 + ?|6|1/2  + 0(6)  , 

(3.5b)  z « A| 6 | 1/2  + 0(6)  , 

where  A are  the  real  roots  of  (2.22b).  They  are  given  by, 


(3.6)  A - A1 (5)  = 


+80 (sgn6) 
40 


Hence  there  are  two  corresponding  inner  expansions  z~  for 
6 > 0,  because  A4  are  real  for  all  values  of  £.  They  are 
shown  by  dashed  curves  in  Figure  2a. 

If  6 < 0,  then  for  real  A4,  5 is  restricted  to  the  inter- 
vals |5|>?0  = 2/T£.  Thus  for  6 < 0,  there  are  four  inner  ex- 
pansions, A4  for  £ > 5°,  and  for  £ < £°.  They  are  shown  by 
dashed  curves  z4  in  Figure  2c.  At  5 “ S0*  both  z+  and  z~  join 
and  have  vertical  tangents,  thus  forming  a "nose”  of  these 
perturbed  branches. 

At  the  limit  point,  we  have  a ■ b ■ 2 and  from  (2.25b) 


(3.7) 


X = XA  + Cl « | + 0 (62) 
z = -B  + A(C)  I « 1 1/2  + 0(6)  , 


where,  the  roots  A of  (2.22a)  are  given  by, 

(3.8)  A = A1  * ±I^-(sgn6)XJl/6]1/2  . 

1 

They  are  real  if  and  only  if  £ £ = -g-  (sgn6).  Hence  there 

are  two  corresponding  inner  expansions,  z~.  They  are  shown 
by  dashed  curves  in  Figures  2a, c.  From  these  figures,  we  ob- 
serve that  as  6 increases  (decreases)  from  zero,  the  limit 
point  moves  to  the  right  (left) . 

For  this  problem  it  is  obvious,  and  it  is  shown  in 
Figures  2a  and  2c,  how  the  inner  and  outer  expansions  connect. 
We  will  therefore  merely  verify  that  the  matching  conditions 
(2.32)  are  satisfied.  We  consider  only  the  bifurcation  point. 
A similar  analysis  applies  to  the  limit  point.  We  insert 
(3.5a)  and  6 = (sgn6)y  into  the  outer  expansions,  and  expand 
the  result  in  a power  series  in  y,  to  obtain 


y*0)  - S2ip-  w + 0(y2)  * nxy  + 0(y2) 


(3.9) 


' »U>  “ (“f  + A)U  + 0(v2)  - "IK  + Of2)  - 


where  we  have  used  the  notation  of  (2.35).  The  inner  expan- 
sions (3.5),  and  (3.6)  at  the  bifurcation  point  give,  for 


large  (, 


* 
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z+  ■ -h + ir*  f + + + 0(^2) 

(3.10) 

*“-^[-^+--:>  + 0(v2>  • 

The  matching  condition  (2.37)  with  k = 1 states  that 

(3.11)  lim  (n.-z.)  = 0 . 

m~  1 1 

Therefore,  it  follows  from  (3.9)- (3.11)  that  z+  must  match 
with  y^j  as  K m (X  ®) , and  it  must  match  with  y*^j  as 
5 (X  -•») , as  shown  in  Figures  2a  and  2c.  Similarly, 

z~  must  match  with  y*Qj  and  y*“j  as  ( + »,  £ -*■  respectively. 

The  composite  expansions  formed  from  the  inner  and  outer 
expansions  then  give  the  perturbed  branches,  as  shown  in 
Figures  3.  This  summarizes  the  analysis  of  Equation  (3.1). 

The  curves  in  these  figures  represent  cuts  6 *»  constant  in  the 
solution  surfaces  of  (3.1)  in  y,  X,  6 space. 

The  Canonical  Cubic. 

For  6 « 0,  the  limit  point  and  the  bifurcation  point 
coalesce  to  form  a new  singular  point  at  X = 1,  y * 0.  Then 
(3.1)  is  reduced  to, 

(3.12)  F[y;X,6]  - y3  + (l-X)y  - X6  - 0 . 

The  bifurcation  branches  are 

(3.13)  y - Yq(X)  5 0,  y Y*(X)  - ±/T=T  . 
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They  are  sketched  in  Figure  4a.  Thus  there  are  four  bifur- 
cation sub-branches.  Wte  observe  that  if  y(X,6)  is  a solution 
of  (4.1),  then  -y(X,-6)  is  also  a solution.  Thus  we  need 
only  analyze  (4.1)  with  6 > 0. 

The  four  outer  expansions  are 

*(0)  ' rr  { ♦ 0({2> 

(3.14) 

y*1}  - ±/T=I  " « + 0 ( 6 2 ) . 

They  are  unbounded  at  the  bifurcation  point,  as  we  show  in 

Figure  4b.  To  obtain  the  inner  expansions,  we  observe  first 

that  F “ 6z  and  F „ - 6.  Thus  from  (2.23)  and  (2.30), 
zz  zzz 

Q « 0 and  C j*  0 and  the  nonlinearity  is  cubic  at  the  bi- 

furcation point.  Furthermore  F°  *=  0,  so  that  A ■ 0.  Thus 
from  (2.31)  and  (2.33)  we  have, 

(3.15)  6 (y)  » u3  , X » 1 + C|6|2/3  + 0(6) 

(3.16)  z - A(C) |«|1/3  + 0(62/3) 

where  the  normalized  eigenfunction  is  ^ « 1,  and  from  (2.29b), 
A is  a root  of  the  cubic 

(3.17)  A3-FA-1«0. 

Our  method  yields  only  a minor  simplification,  since  (3.12) 
is  essentially  a canonical  form  of  the  algebraic  equation 
(2.29b)  to  which  general  problems  with  cubic  nonlinearities 
are  reduced. 
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It  is  easy  to  show,  by  using  elementary  properties  of 
cubic  equations,  that  (3.17)  has  only  one  real  root 

(3.18)  A+(£)  > 0 , if  £ < £°  = 3/  (2) 2/3  , 

and  it  has  three  real  roots 

(3.19)  A+(£)  > 0 , a"1}(£)  < 0 » AU)(t)  < 0 ' 

if  £ > 5°  . 

Thus  there  is  one  inner  expansion  for  £ < £°  and  three  inner 
expansions  for  £ > £°.  The  matching  conditions  (2.37)  decide 
how  the  inner  and  outer  expansions  connect. The 
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results  of  the  matching  are  summarized  in  Table  I. 
Table  I.  Inner  and  Outer  Matching  as  6 -»  0. 


The  resulting  inner  expansions  z+,  z^j , z"j2j  are  shown  by 
the  dashed  curves  in  Figure  4b.  The  composite  expansions  of 
the  perturbed  branches  are  shown  in  Figure  4c. 

A similar  analysis  can  be  applied  to  the  cubic 

(3.20)  F(y;X,6]  « y3  - (l-X)y  - X6  - 0 . 

The  resulting  bifurcation  branches,  inner  and  outer  expan- 
sions, and  the  composite  expansions  are  shown  in  Figures  5. 

We  observe  that  there  are  bifurcation  branches  for  X < X ■ 1. 
This  subcrltical  bifurcation  is  in  contrast  to  the 

supercritical  bifurcation  exhibited  by  (3.12)  with  6 ■ 0, 
where  the  bifurcation  branches  exist  only  for  X > XQ.  The 
bifurcation  from  X ■ 1 of  the  solutions  of  (3.2)  is  another 
example  of  subcrltical  bifurcation.  The  effect  of  small  in- 
homogeneities is  most  dramatic  for  subcrltical  bifurcation  , 
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4.  A Nonlinear  Boundary  Value  Problem. 

We  now  apply  the  Method  to  the  nonlinear  boundary  value 
problem, 

(4.1)  y”  ♦ X (G(y)+4g(x,y)  ] ■ 0 , y{0)  - y(w)  * 0 , 

for  the  function  y(x).  Here,  G(y)  and  g(x,y)  are  prescribed 
functions.  Furthermore,  we  assume  that  they  satisfy  the  con- 
ditions 

(4.2a)  G (0)  - 0 , 6*(0)  s ax  > 0 f g(x,0)  f 0 , 

and  that  G(y)  has  the  asymptotic  expansion, 


The  bracketed  term  in  (4.1)  can  be  interpreted  as  a source, 
or  growth  term,  where  X is  proportional  to  the  strength  of  the 

source.  For  4 ■ 0,  the  source  is  G(y).  Then  4g  can  be  con- 

I 

sidered  as  an  "impurity”  in  the  source,  where  4 characterizes 
the  magnitude  of  the  impurity.  If- 
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g is  independent  of  y,  then  fig(x)  is  equivalent  to  an  ap- 
plied external  forcing  function.  Since  the  sign  of  5 can  be 
absorbed  in  g , we  shall  assume,  with  no  loss  of  generality, 
that  6^0. 

The  bifurcation  problem  corresponding  to  (4.1)  is 

(4.3)  y"  + XG(y)  « 0 , y (0)  » y(ir)  - 0 . 

Since  G(0)  » 0,  y - Yq  = 0 is  a solution  of  (4.3)  for  all 
values  of  X.  The  singular  points  of  this  solution  branch, 
which  are  all  bifurcation  points,  are  the  eigenvalues  XR  of 
the  linearized  problem  (2.2)  corresponding  to  (4.3), 

(4.4)  ♦"  + Xa^  - 0 , f(0)  «=  *(ir)  = 0 . 

They  are  given  by 

(4.5a)  Xn  ■ n^/a^  , n ■ 1,2,...  • 

The  corresponding  eigenfunctions  are 

(4.5b)  4n  ■ /27»  sin  nx  , n ■ 1,2,... 

We  determine  the  solution  branches  of  (4.3)  that  bifur- 
cate from  X^  by  a perturbation  expansion.  This  expansion  is 
valid  for  X near  X^.  Similar  expansions  can  be  obtained  for 
the  branches  bifurcating  from  Xj , X^ , • • • . Thus  the  bifurca- 
tion branches  are  known  only  near  the  bifurcation  points,  and 
then  only  approximately.  The  analysis  of  (4.1)  will  thus 
demonstrate  the  application  of  our  method  to  such  problems. 


b 
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To  obtain  the  perturbation  expansions  for  the  bifurca- 
tion  branches,  we  define  an  amplitude  parameter  e. 


(4.6) 


€ = <y,4,> 


where  the  inner  product  is  defined  by, 

IT 


<f 


,g>  = | f (x)g (x)dx  . 
0 


We  seek  expansions  of  the  solutions  of  (4.3)  near  y = 0, 


X « X^  in  the  form 


(4.7)  y(xje) 


Yj  (x)  , X ® Xj  + . 


Then  by  substituting  (4.7)  into  (4.3)  and  (4.6),  wo  obtain, 
in  the  usual  way,  a sequence  of  linear  boundary  value  prob- 
lems to  determine  the  coefficients  in  (4.7).  An  analysis  of 
these  problems  determines  the  leading  terms  in  these  expan- 
sions as 


(4.8a) 

(4.8b) 

(4.8c) 


y(x;e)  » e^fx)  + 0(0 

a 
*1 


X ■ X,  - 


* jfr  7 + 0<e‘>  ’ if  “2  * 0 


X - X,  - 


e2  — + 0(c3)  , if  a,  ■ 0 and  a,  / 0 
2ira^  J 


We  denote  the  bifurcation  sub-branches  that  are  approximated 
by  (4.8)  by  y ■ Y*(xje).  We  observe  that  the  bifurcation  is 


i __ 
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subcritical  if  a2  + 0 and  is  supercritical  (subcritical)  if 
a2  ■ 0 and  a3  < 0 (>  0) . 

The  Perturbed  Problem  (&  y*  0) . 

The  outer  expansions  of  the  perturbed  problem  are  ob- 
tained by  inserting 

m 

(4.9)  y(x;X,«)  - Yj  (x;  1)  6^ 

into  (4.1),  where  yQ(x»X)  is  any  solution  sub-branch  of  the 
bifurcation  problem  (4.3).  This  gives 

(4.10a)  F°yjL  = y[  + XG^y^yj^  - -Xg(x,yQ)  s RjtxJX) 

yx(0;X)  - y1(ir?X)  « 0 , 

(4.10b)  F°y.j  - Rj  (x; X)  , Yj(0;X)  - Yj(ir;X)  - 0 , 

j * 2,3,0..  m 

Here  the  inhomogeneous  terms  Rj  are  determined  by  Yq'Yj #•••* 

Vj-I- 

We  solve  these  linear,  boundary  value  problems  by  using 
eigenfunction  expansions.  We  employ  the  eigenfunctions  vm 
of  the  operator  F°.  That  is,  the  v^  are  determined  by 

(4.11)  v‘  ♦ >G’(y0)vm  - 8nvn  - 0 , vnm  - Vm(w)  - 0 , 

where  a are  the  corresponding  eigenvalues.  Then  by  expand- 
m 

ing  the  inhomogeneous  terms  Rj  as 
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(4.12) 


we  obtain  the  solutions  of  (4.10)  as 


(4.13)  Yj  (x»  X)  *= 


Jjjn 


(M 


o„UT 


v ( x ; X ) , j a 1 » 2 , . . . 


m 


m 


A ^ 


Por  the  branches  yQ  «*  Y = 0 , we  have  G ' (yQ ) = a^ , and 
the  eigenfunctions  and  eigenvalues  of  (4.11)  are. 


(4.14)  vm  - 4m  , om  - Xax  - in2  - a^X-XJ  , 

in®  If  2|  • • • • 

We  denote  the  outer  expansions  (4.9)  corresponding  to  Y*  by 
y*0j.  Then  by  using  (4.13)  and  (4.14)  we  obtain, 

m 

(4.15)  y*0) . , 

where  gm  is  defined  by 

X 

(4.16)  gm  5 | g(x,0)4m(x)dx  , m - 1,2,...  . 

0 

In  the  notation  of  Section  2,  g^  is  proportional  to  I,  and 

hence  by  assumption  (2.8),  g^  f 0.  Thus  as  X -*■  X^  from  below, 

y*Qj  is  positive  (negative)  if  g1  > 0 (<  0).  The  signs  of  y*Qj 

are  switched  if  X ♦ X^  from  above.  The  perturbed  branches 

(4.15)  are  unbounded  at  the  bifurcation  point. 

For  the  bifurcation  branches  y^  * Y^(x;e),  X ■ X(e), 

given  by  (4.8),  the  coefficients  of  the  differential  equations 


I 
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in  (4.10a)  and  (4.11)  depend  on  x and  the  parameter  e.  Thus 
we  solve  the  eigenvalue  problem  (4.11)  by  a perturbation  ex- 
pansion in  e.  That  is,  we  seek  solutions  of  (4.11)  in  the 
form 


(4.17)  vm(x,e)  - >.  vmj(x)e oje) 


°m,je  ' 


where  <vjn>  ■ 1.  An  analysis  of  the  resulting  system  of  linear 

boundary  value  problems  for  the  v . yields  the  following: 

m,j 

(4.18)  Vm“*m+0<E2>  t m-1,2,...  ; 


for  a9  / 0, 


(4.19a)  o. 


X - 4(1) 3/2  ^ 


) , om  - (1-m  ) + 0(e)  , 


in  * 2,3,.*.  I 


for  ■ 0,  and  a~  f 0, 


(4.19b)  ax  - e2  + 0(e3)  , om  - (1-m2)  + 0(e3)  , 

m ® 2,3,...  . 

To  evaluate  the  coefficients  in  the  outer  expansions  cor- 
responding to  the  bifurcation  branches  from  (4.13),  we  require 
the  Fourier  coefficients  m.  For  example, 

(4.20)  Rj^  m ■ | * "*<*)  | g(x,Y^  (x;e)  Jvm(x;e)dx  . 

He  insert  the  expansions  (4.8),  and  (4.18)  into  (4.20),  to  obtain 
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(4.21) 


.-1 


1,«  " -1  9m  * 0(c) 


Finally,  by  substituting  (4.8),  (4.18),  (4.19)  and  (4.21) 
into  (4.13),  and  then  substituting  the  result  into  (4.9),  we 
obtain  the  outer  expansions,  which  are  denoted  by  y^,  as 


(1) 


2 f3  r-n- 3 3/2  g.  -i 

e^x)  ♦ 0 (e  j - [f|JJ  5-F  *1  + °<1)J6 


+ 0(«2)  , 

if 

a2  / 0 

(4.23) 

y*1(  - «*XW  + 0U2)  - [(|) 

a3e 

_ A,  A 

0 (l/cTj 

r t ' 

+ 0(«2)  , 

if 

a2  * 0 

and 

a3  * 0 

We  observe  that  the  perturbed  branches  (4.23)  are  unbounded 
at  the  bifurcation  point  X^  (e  * 0) . 


The  Inner  Expansions. 

The  second  functional  derivative  of  the 
nonlinear  operator  (4.1)  at  the  singular  point  X ■ X^,  y ■ 0 
is  f®2  ■ 2X^a2>  Therefore,  Q given  by  (2.23)  is  propor- 
tional to  a2  , and  the  boundary  value  problem  (4.1)  has  a 
quadratic  nonlinearity  if  a2  f 0.  If  a2  ■ 0,  and 

a^  ¥ 0,  then  it  is  easy  to  see  that  C given  by  (2.30)  is  not 
sero  and  the  boundary  value  problem  has  a cubic  non- 

linearity at  the  bifurcation  point. 

First,  we  shall  study  (4.1)  with  a2  f 0.  Since  F°  = 0, 
and  hence  A - 0,  the  analysis  of  Section  2 for  6 > 0,  shows 


-r 


that  we  must  seek  solutions  of  the  perturbed  problem  in  the 
form 


(4.24)  z (x;  i0  - y (xj  X(y)  ,6  (p) ) ° , z..  (x)u3  , 


and  X « X.  + 54PJ  # 6 » p . 


The  coefficients  z^  then  satisfy, 

(4.25)  + Zj  * px  s 0 , zx(0)  - z^z)  « 0 , 

o 2 

(4.26)  Fzz2  - p2  H - gU,!))^  , 

0 

z2(0)  * z2(z)  «=  0 , 


(4.27) 


rlzi  - “j  • 


Zj  (0)  *=  Zj  (z)  -0  . 


Thus  Zy , is  given  by 

(4.28)  z^  - A4»1(x)  , 

where  A satisfies  the  quadratic  equation 

(4.29)  A2  + 2p(A  + q - 0 , 


3fz)2^2  al 

where  2p  = |[7]  — , q 


3fz)3/2  gl 

IT  aT 


There  are  two  real  roots  of  (4.29)  if, 


(4.31)  C2  > (q/p2)  . 
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Thus  there  are  two  inner  expansions 

(4.32)  z*  * A* (£) sin  x 61/2  + 0(62)  , 

X - \x  + £61/2  +0(6)  . 

If  *2gi  < then  these  expansions  are  defined  for  all  £.  If 
a29l  > t*ien  there  are  two  expansions  for  £ > q*^2/p,  and  two 
expansions  for  C < - q 1^2/p.  At  |£|  « q1/,2/p,  IdA*/^)  * «. 


The  Matching. 

To  apply  the  matching  condition,  we  first  express  the 
outer  expansions  in  terms  of  the  inner  variables  as  in  (2.35). 
Thus  substituting  (4.24)  into  (4.13),  (4.15),  and  (4.16),  and 
(4.23),  and  rearranging  the  result  in  a power  series  in  p,  we 
get. 


(0)  V0'-2’ 


(4.33) 


- -[«  - rl]V  + 0 (p2)  + ...  . 
The  asymptotic  forms of  the  inner  coefficients  are. 


(4.34) 


lim  zt  - lim  *7  - “ (x)  , 


lim  z7  - lim  z+  m -2pW,  (x)  . 

E-w»  A F+-m  1 1 


t 
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Since  is  the  coefficient  of  y in  (4.33),  it  follows  from 
(4.33)  and  (4.34)  and  the  matching  condition  (2.37)  with 
k • 1 that  the  branches  connect  as  summarized  in  Table  II. 
The  resulting  composite  expansions  for  g^^  > 0 and  g1  < 0 are 
shown  in  Figures  6. 


Table  II 


h ♦ 

00 

—00 

+ 

+ 

z 

o 

y(l) 

+ 

z 

y(l) 

o 

>T 

The  analysis  of  (4.1)  with  a2  *»  0 and  a^  j-  0 (the  cubic 
nonlinearity)  is  similar,  and  hence  we  present  the  results 
only.  The  stretching  variables  are  now 

(4.35)  X - \x  + Kv2  + 0 (y4)  , 6 - y3  . 

The  inner  expansions  are,  c.f.  (2.27) 

z - A4>jy  + 0 (y2)  , 

where  A is  any  real  root  of  the  cubic,  c.f.  (2.29b) 

(4.36)  A3  ♦ 2r£A  ♦ 8 - 0 , 

2 IT 

where  ^ H ■ 
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Thus  if  dj  < 0 (r  < 0) , which  corresponds  to  supercritical 
bifurcation,  we  obtain  essentially,  the  cubic  (3.17).  The 
graphs  of  the  resulting  composite  expansions  are  similar  to 
Figure  4c.  If  a^  > 0 (r  > 0) , so  that  we  have  subcritical 
bifurcation,  the  analysis  is  equivalent  essentially  to  that 
of  (3.22).  Thus  the  graphs  of  the  resulting  composite  ex- 
pansions are  qualitatively  similar  to  Figure  5c. 

Interpretation . 

To  interpret  the  solutions  of  the  perturbed  problem,  we 
consider  (4.1)  as  describing  the  equilibrium  states  of  a 
diffusion  process,  so  that  we  study  only  non-negative 
solutions.  The  stability  of  these  equilibrium  states  is 
determined  from  the  linearized,  time  dependent,  diffusion 
problem  that  is  linearized  about  a solution  y(x)  of  (4.1). 

First  we  consider  (4.1)  with  a2  > 0.  The  composite  so- 
lutions are  shown  in  Figures  6.  For  the  bifurcation  problem 
(6  “ 0),  the  zero-concentration  solution,  y = 0,  is  stable 
for  X < and  unstable  for  A > A^.  Thus  for  the  bifurcation 
problem,  as  X (the  source  strength) 

exceeds  X^,  the  zero-concentration  state  is  unstable,  and  a 

transient  diffusion  occurs.  For  the  perturbed  problem  (6  0) 

with  g^  > 0,  Figure  6a,  as  X increases  from  zero  there  is  a 
"small-concentration"  equilibrium  state.  When  X exceeds  Au, 


a 
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Appendix  I.  A Model  of  Elastic  Arch  Buckling. 

Two  rigid  and  massless  rods  of  length  t are  connected 
by  a linear  torsional  spring,  as  shown  in  the  Figure  7. 

The  total  mass,  m,  of  the  system  is  concentrated  at  the  junc- 
tion point.  The  left  and  right  endpoints  of  the  rods  are 
pinned,  so  that  they  cannot  move  vertically,  but  they  are 
free  to  slide  horizontally.  A compressive  force  T is  applied 
to  the  ends  of  the  rods,  as  shown  in  the  figure.  The  initial 
position  of  the  rods,  before  application  of  the  force  is 
shown  by  the  dotted  lines.  Thus  in  the  initial  position,  the 
junction  point  is  a specified  distance  wo  above  the  horizontal 
and  the  rods  have  an  initial  slope  0Q.  As  a result  of  the 
force  the  junction  point  is  a distance  w above  the  horizontal 
and  the  deformed  slope  is  6. 

The  motion  of  the  rods  is  resisted  by  the  torsional 
spring,  which  has  a spring  constant  T,  and  a nonlinear  exten- 
sional  spring.  Then  the  torsional  resistance  moment  M,  is 
proportional  to  the  change  in  angle  of  intersection  of  the 
rods.  Thus  we  have, 

m - ne-eo)  . 

The  force  in  the  extensional  spring  is  a function  of  the  re- 
lative displacement  w - wQ.  In  particular,  we  shall  consider 
the  nonlinear  law 

2 3 

F * a(w-w  ) + b(w-w  ) + c(w-w,J  , 

o o o 


% 


where  a,  b and  c are  the  prescribed  spring  constants.  If 
b ^ 0,  then  the  spring  is  unsymmetrical. 

Since  sin  0 * w/A  and  sin  0O  * wq/A,  we  have  for  small 


0 and  0Q, 


0 = w/A  , 


eo  : »e/i  . 


Then  by  defining  a dimensionless  variable  y by 
y = ^2  , A 5 l(a+4r/A2)/c]1/2  , 


the  equation  of  motion  of  the  mass  is, 
(A.l)  y + y3  + 2{$y2  + (l-X)y  - XS  « 0 


The  variable  y is  a dimensionless  "excess"  deflection  of  the 
mass.  The  constants  6,  X and  6 in  (4.1)  are  defined  by 


(A. 2) 


« = b 
6 - 7S5  ' 


x _ 2T 


X • 


The  parameter  X,  which  is  proportional  to  the  thrust,  is 

called  the  load  parameter,  and  6 is  the  imperfection  para- 

— 

meter.  For  equilibrium  states  y ■ 0,  and  Equation  (A.l)  is 
reduced  to  (3.1). 


, 
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Captions  for  Figures 


Figure  1.  None 

Figure  2a.  The  inner  and  outer  solutions  for  the  perturbed 
Equation  (3.1),  with  6 > 0.  The  inner  branches  are  repre- 
sented by  the  dashed  curves,  and  the  outer  branches  by  the 
heavy  solid  curves.  The  light  solid  curves  represent  the 
bifurcation  branches.  In  these  figures,  we  have  taken  B in 
the  interval,  0 < 0 < 1. 


Figure  2b.  The  bifurcation  branches  (6  ■ 0)  for  the  elastic 
arch  model.  Equation  (3.1).  Y*  are  the  sub-branches  of  YQ, 
and  and  Y^~  are  the  sub-branches  of  Y^. 


Figure  2c.  The  inner  and  outer  branches  for  Equation  (3.1) 
with  6 < 0. 


Figure  3a.  The  composite  solutions  for  Equation  (3.1)  with 

6 > 0. 

Figure  3b.  The  composite  solutions  for  Equation  (3.1)  with 
6 < 0 and  C°  ■ 2/Jff. 


Figure  4a.  The  bifurcation  branches  (0  • 0)  for  Equation 
(3.12). 


i 


Figure  4b « The  inner  and  outer  solutions  for  Equation  (3.12) 


with  6 > 0. 

Figure  4c.  The  composite  solutions  for  Equation  (3.12)  with 
6 > 0.  The  graphs  of  the  composite  solutions  for  Equation 
(3.12)  with  £ < 0 are  the  images  with  respect  to  the  X axis 
of  the  curves  shown  in  the  figure. 

Figure  5a.  The  bifurcation  branches  (6  «•  0)  for  Equation 
(3.22). 

Figure  5b.  The  inner  and  other  solutions  for  Equation  (3.22) 
with  6 > 0. 

Figure  5c.  The  composite  solutions  for  Equation  (3.22)  with 
£ > 0.  The  images  of  these  curves  with  respect  to  the  X 
axis  are  the  graphs  of  the  composite  solutions  for  (3.22) 
with  £ < 0. 

Figure  6.  The  composite  solutions  for  the  nonlinear  boundary 
value  problem  (4.1)  with  a2  > 0.  The  noses  of  the  branches 
for  gx  > 0 are  Xu  - Xx  - (q^2/p)6^2  + •••,  and 
XL  - Xx  ♦ (q1/2/P>  fi1/2  *■  • • • . 

Figure  7.  Model  of  an  elastic  arch. 


